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Abstract 

A continuum-QRPA approach to calculation of the 2v/3/3- and Oi//3/3-amplitudes has been for- 
mulated. For 130 Te a regular suppression (about 20%) of the high-multipole contributions to the 
Oz^/3/3-amplitude has been found which can be associated with additional ground state correlations 
appearing from the transitions to collective states in the continuum. At the same time the total 
calculated Oi//3/3-amplitude for 130 Te gets suppressed by about 20% as compared to the result of 
the usual, discretized, QRPA. 

1 Introduction 

At present, the most elaborated analysis of the uncertainties in the Ozz/^-decay nuclear matrix elements 
M 0u calculated within the QRPA and RQRPA has been performed in the recent works QjJ Ej. The 
experimental 2vf3~ f3~ -decay rate has been used to adjust the most relevant parameter, the strength 
g pp of the particle-particle interaction. With such a procedure the values of M 0u have been shown to 
become essentially independent on the size of the single-particle (s.p.) basis. Furthermore, the M 0u, s 
have been demonstrated to be also rather stable with respect to the possible quenching of the axial 
vector strength as well as to the uncertainties of parameters describing the short range nucleon 
correlations. 

The calculations in E] have been performed within "the standard QRPA" scheme in which the 
BCS ground state and the spectrum of the excited states are built on a discrete s.p. basis. However, 
in order to really assess the uncertainty of the QRPA calculations of the M 0v one has to address 
the questions regarding the accuracy of "the standard QRPA" itself, and whether there are not some 
important contributions missing. 

The accepted procedure [TJ E] of averaging the results of the calculations performed with different 
choices of the model space size looks rather weakly justified. One can expect a priori that enlargement 
of the model space should lead to more accurate results (in other words, any basis truncation leads to 
an uncertainty). The usual statement that "in the QRPA one can include essentially an unlimited set 
of single-particle states" [T] can be considered at presence only as principally correct, but not realized in 
practice. In reality, the large number of the major shells ^> 1 in the QRPA calculations can only be 
achieved by adding low-lying major shells composed of the bound states. Basically, only one major shell 
lying higher than the Fermi-shell one can be considered, because one immediately encounters principal 
limitations of the approximation of the single-particle continuum by discrete levels. Neglecting the 
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single-particle continuum leads to a missing strength in the usual QRPA calculations, especially for 
description of the high-multipole excitations with L > 2 1 . 

The contribution of these multipoles to M 0u is particularly important because the monopole (Fermi + 
Gamow- Teller) one is suppressed due to the symmetry constrains (see, e.g., the multipole decomposition 
of M 0u in Fig. 5 of |2]; for a recent discussion how the SU(4)-symmetry violation by the particle-particle 
interaction affects M 2v see |3]). Furthermore, the pure pairing contribution to M 0u is almost completely 
(by an order of magnitude) suppressed by the ground state correlations, short-range correlations etc., 
therefore fine effects can be expected to play an important role. 

The question about the dependence of the QRPA results on the s. p. -basis size as a source of the 
uncertainties in the calculated M 0l/ 's would be completely resolved if one could include the entire s.p. 
basis into the calculation scheme. The only possible way principally to perform such ultimate-basis 
QRPA calculations is provided within the continuum-QRPA. Also, to have an alternative formulation 
of the QRPA can help to promote our understanding of the current QRPA results and their deficiencies 
to a higher level. In particular, the continuum-QRPA provides a regular way of using realistic wave 
functions of the continuum states in terms of the Green's functions and there is no need to approximate 
them by the oscillator ones. 

Two principal effects of the inclusion of the s.p. continuum within the pn-QRPA, which affect M 0u, s 
in an opposite way, can be expected. First, pairing in the continuum can increase respective 0^/3/3 
sum rules (the smallness of the A/E in the continuum can be compensated by a large corresponding 
partial s.p. matrix element). Second, additional g.s. correlations can appear due to collective multipole 
states in the continuum that decreases M 0v . It is noteworthy that at the moment the continuum-QRPA 
consistently including pairing in the continuum has not been formulated yet and only the continuum- 
QRPA with the pairing realized on a discrete basis can be used. 

The principal aim of this work is to formulate for the first time a method to calculate the double 
beta decay matrix elements within the continuum-QRPA. 

2 Continuum-QRPA 

The continuum- RPA has been used for a long time to successfully describe structure and decay properties 
of various giant resonances and their high-lying overtones embedded in the single-particle continuum. 
However, to apply the continuum-RPA in open-shell nuclei one has to take the nucleon pairing into con- 
sideration and develop a continuum-QRPA approach. The approach should account for the important 
effects of particle-particle interaction along with the usual particle-hole one. Such a continuum-QRPA 
approach has been developed in |U Ej. In |3] the approach has been applied to the analysis of the 
low-energy part of the Gamow- Teller (GT) strength distribution for description of the single-beta decay 
relevant to astrophysical applications. 

2.1 QRPA equations in the coordinate representation 

The system of homogeneous equations for the forward and backward amplitudes X^J^ and Y^ s \ 
respectively, is usually solved to calculate the energies uj s and the wave functions | J/i, s) of the isobaric 
nucleus within the pn-QRPA (see, e.g., [B]). It is impossible to handle the infinite number of the 
amplitudes X, Y if one wants to include the single-particle continuum. Instead, by going into the 
coordinate representation the pn-QRPA can be reformulated in equivalent terms of four-component 
radial transition density {o\ J (I = 1, ... 4). The components are determined by the standard QRPA 
amplitudes X and Y as follows: 

1 Note, that the QRPA is barely suitable for the description of the multipole contributions to M 0l/ with L > 5 because 
they are completely dominated by the short-range behavior of the wave function. 
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where u, v are the coefficients of Bogolyubov transformation, x* 



*(3m Xn{r)xAr) with t$ m 



^j^ ('k\\Tjls\\ 1/ ) being the reduced matrix element of the spin-angular tensor Tjls^ and x*( T ) {Xv( r )) 
being the radial wave function of a single-particle proton (neutron) state. Hereafter we shall systemat- 
ically omit the superscript "J^s" when it does not lead to a confusion. 

According to the definition (JTJ), the elements Qi, Q2, Q3, 84 can be called the particle-hole, hole- 
particle, hole-hole and particle-particle components of the transition density, respectively. In particular, 
the transition matrix element to the state \s, Jfi) corresponding to a probing particle-hole operator 

» ' \ T/"_ ( rr> \ ' I 1 „ n~ 1 1 1 r\ r\^~ nr l m r 4" Vi i"i «">1 r>m r~\t~\\- /~\ . • f /^'^ ' ^ / 
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r a is determined by the element g x : J q\ ' (r)Vj(r) dr. 



is as follows: 



The pn-QRPA system of equations for the elements g\ 

' A i I J p(rr',u; = uj s ) F { p(r'r") g ( f s \r") dr'dr", 



q\ >(r) 



E 

K 



or schematically, denoting all the integrations and summations as {. . .}, g — {AFg}, where '(r^j 
represents the residual interaction in .^-channel (K = 1,2 — p-h channel, K = 3, 4 — p-p channel) 
after separation of the spin-angular variables. The 4x4 matrix Aix{ r i r 2i w) is the radial part of the free 
two-quasiparticle propagator. The expressions for the elements of the free two-quasiparticle propagator 
Ajk can be obtained by making use of the regular and anomalous single-particle Green's functions for 
Fermi-systems with nucleon pairing in an analogous way to how it was done in the monograph [7] to 
describe the Fermi-system response to a single-particle probing operator acting in the neutral channel |5J. 
The corresponding analytical representations for the elements v4/^(r 1 r 2 , co»)are: 



A IK (nr 2 ,u) = e x,An)x,u(r 2 ) AT K {w)\ Ay K = A^ (3) 

TTU 

_ ujvl ujvl _ u\ul v 2 v vl 

_ , . TP TP , . TP TP ' ^33 — , , 77 rp + 



uj — E n — E v — uj — E n — E v uj — E n — E v — uj — E n — E v 

12 uj — E n — E v —uj — E^ — Ey 

Aw - 7/ 7 , l U ^ V * ) A*" - o, v I ^ 

A VA — u v v v\ „ .tp TP >■> A U ~ UnV7T 



uj — E n — E v — uj — E n — E v uj — E n — E v — uj — E n — E v 

A™{uj) = A™{-uj), AH(u) = A%(-u>), A%(u>) = A%(-u>), A%(u>) = A%(-u>). 

The total two-quasiparticle propagator A that includes the QRPA iterations of the p-h and p-p 
interactions satisfies a Bethe-Goldstone-type integral equation: 

A = A + {AFA} (4) 

and has the following spectral decomposition: 

x&\nw) = e gr " (r ' )g *"' w - E ( 5 ) 

g UJ — UJ S + ZO g UJ + UJg — ZO 

Thus, all the necessary information about the QRPA solutions resides in the poles of A^ J7T \ 
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2.2 Strength functions 

One can use Im A to calculate different strength functions. Strength function corresponding to a charge- 
exchange single-particle probing operator 

Vff = Y.y^a)Tj L s,{n a )T^ (6) 

a 

acting in /3^-cnannel is defined by the usual expression: 

5 (T) M = E|(^IK7 ( ?|o)| 2 ^-^) (7) 

s 

with = Ef — E being the excitation energy of the corresponding isobaric nucleus measured from 
the ground state of the parent nucleus. Making use of the spectral decomposition (jSJ) one can easily 
verify the following integral representations of the strength functions in term of Im A 

S^{u~) = -^ImfVj^AiPin^^Vj^dnd^ (8) 

S(+)(o;+) = -HmJVjir^AiPin^^Vj^dnd^ (9) 

or, schematically, 

4" V) = -~ {VA n (u)v} 4 + V + ) = -~ {vA 22 (u)v} 

with uj t = uj ± (A p — A n ). The calculated pn-QRPA spectrum in uj is to be shifted in energy in order to 
be measured from the ground state of the parent nucleus. It has to do with the fact that in the QRPA 
the BCS hamiltonian H — X p Z — X n N is used. 

One can also define a non-diagonal strength function like 

4~ ) («) = e<^ ) i*><^; ) io>%-^) ( io ) 

s 

with uj' s = E s — (Eq + E /)/2. Such a strength function is closely related to the amplitude of the 2vf3[3- 
decay. To calculate S v \uj) within the pn-QRPA one faces the usual problem that the spectrum \s) 
comes out slightly different when calculated starting from the initial or final ground states. Identifying 
the BCS vacuum |0') with that of |0) one gets 

S { -\uj-) = --Im f VjirjAiPinr^^Vj^dndri (11) 

71 J 

or, alternatively, identifying |0) with |0') 

S ( -\io + ) = --Im f Vj^A'^inr^^Vj^dndri (12) 
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where A' is calculated with respect to |0'). 

To calculate the strength functions, it is more convenient to use a system of inhomogeneous cQRPA 
equations in /^"-channel: 

S { y\uj-) = -I Im E / V(n) A 1K ( ri r 2 , uj) V K (r 2 , uj) d n dr 2 , (13) 



K 

7(~ )(,..-\ 



Sf->(u-) = -ilmE JV(r 1 )A 2K (r 1 r 2 ,uj)V K (r 2 ,uj)dr 1 dr 2 , (14) 

K 

V7(r, uj) = V(r)5 n + E / F K (rri) A IK (r x r 2 , uj) V K {r 2 , uj) dr x dr 2 , (15) 

K 



or, in /5 + -channel: 



= - J Im E / V{r) A 2K (nr 2 , uj) V K (r 2 , uj) d n dr 2 , (16) 

K 

Vi(r, uj) = V(r)8 I2 + E J F K (rri) A IK (r x r 2 , uj) V K (r 2 , uj) dr x dr 2 . (17) 

K 
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2.3 Taking the s.p. continuum into consideration 



Up to now the way of taking the s.p. continuum into consideration has not been specified explicitly 
within this formulation of the pn-QRPA. If one lets the double sums in (J3J) run just over the bound 
proton and neutron s.p. states, the version of the pn-QRPA presented in two preceeding sections is fully 
equivalent to the usual discretized one formulated in terms of X and Y amplitudes. We make use of 
this fact and compare discrete-QRPA results calculated in these two different, but formally equivalent, 
ways in order to check the consistency of the scheme. 

To take the s.p. continuum into consideration, one has to do the following in (J3J): 

1. To approximate Vi, Ui and Ei by their no-pairing values v — 0(1), u — 1(0), E — \e — A| for those 
s.p. states which lie far from the chemical potential (i.e. \e — A| ^> A). The accuracy of this 
approximation is t^Au- 

2. To use the s.p. Green function: g( a )( r i r 2, e) = i)Xa r 2) ^ explicitly perform the sum over 
the s.p. states in the continuum. 

As an example of such an approach we present here the final expression for An: 



An(rir 2 u) = V ^ X ^ ri " >X ^ T2 " > 1 12' 12 v lx^(n)x7rA r 2) 

v TV u) - E n - E v v w<7r . to - E v + X n - e n 
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+ 12 {%){*)) Xv(Ti)Xv{r 2 ) vlg[^(r x r 2 , K + w - E v ) 
+ {n ^ u,uj ^ —to} (18) 

where the projected s.p. Green's function is g'^Jr^, e) = 5 f ( 7r )(^i^2, £ ) ~12 — ~~- ^ ne P rrme d 
sum J2' runs over only those s.p. states which comprise the BCS basis (for instance, all proton s.p. 

States 7T m j n < 7T < TTmax)- 

This continuum-QRPA method has been applied in our recent paper [5] to describe the Fermi and 
Gamow- Teller strength distributions in semi-magic nuclei. 



2.4 Description of the (3j3-decay within the cQRPA 

The spectral decomposition of A © can be used for calculation of /3/3-decay matrix elements in a way 
similar to the one described in Sec. 12.21 For instance, the 2^/3/3-amplitude can be calculated according 
to the following expression: 

= ~J A$( ri r 2 ; to = 0) d n dr 2 + 6M% T (19) 
5 MX = ^ f^n^inr^dndr, 



GT n J (u + 5E)u 

where 5E = Q + m e c 2 + X p — X n . We use in deriving (|T9*j) the approximation that the BCS vacuum 
|0') of the final g.s. is taken to be the same as |0) of the initial g.s. The expression (|T9*j) can be further 
rewritten in terms of the effective field Vk (fT5j) : 

m gt = ~ E / 45knr 2 ; to = 0)V K (r 2 ; to = 0) dr L dr 2 + 5M% T (20) 



A' 



7T J (iO + 8E)UJ 



The same procedure can be applied to calculate within the cQRPA the matrix element of a two-body 
operator 

Vt~ ] = E E Vj L (r a , r 6 )T JL5M (n a )r; iSM (n 6 )ri-)r 6 ( - ) (21) 

ab JLSfj, 

between the ground states |0) and |0') as a sum of all partial contributions M^ JV) : 

MS— ) = (0'\vt~ } \0) = ^M (Ji) (22) 

JL 

M (j L ) = _(m±lj duJ fVjLir^^ImAPirx^^dnd^ 

(the identification of the ground states described above has to be done). 

The neutrino potential V% "* in the simpliest (but rather rough) Coulomb approximation has the 

well-known partial radial components VjL(ri,r 2 ) = ^Tir^T (j^J ( r < = m i n { r li r 2)-, r > = max(ri,r 2 )). 
When the Jastrow factor (to account for the short range correlations) and the energy dependence of 
the neutrino propagator are considered, the decomposition of the neutrino potential over the Legandre 
polinomials (|2*T|) can be done numerically. 

3 First calculation results 

For the first calculations of M 2u and M 0u within the continuum-QRPA we adopt a rather simple 
nuclear Hamiltonian similar to that used in [HUH]- The chosen nuclear mean field U(x) consists of the 
phenomenological isoscalar part U (x) along with the isovector U\(x) and the Coulomb Uc(x) parts, both 
calculated consistently in the Hartree approximation (see The residual particle-hole interaction as 
well as the particle-particle interaction in both the neutral (pairing) and charge-exchange channels are 
chosen in the form of the zero-range, 5-functional, forces (hereafter all the strength parameters of the 
residual interactions are given in units of 300 MeV-fm 3 ). 



Table 1: The fitted parameters /i, oL and the calculated M 0u for 130 Te (g^ = 1.25). 



QRPA 


fph 


9pp 


M 0u 


discrete 
continuum 


0.60 
0.65 


1.20 
1.15 


2.24 
1.79 



Fixing the model parameters is done as follows: 

• The pairing strengths g^ atT , g^ atT are fixed within the BCS model to reproduce the experimental 
pairing energies. 

• The p-h isovector strength f® h is chosen equal to unity, f® h = 1.0 that allows to reproduce the 
experimental nucleon binding energies for closed-shell nuclei provided the isospin-selfconsistency 
of the isovector p-h interaction and the symmetry potential U\(x) of the mean field is used (see |5j). 

• The p-h spin-isovector strength f£ h is fitted to reproduce the experimental energy of the GTR. 

• By choosing the p-p isovector strength <j£ = (gP mr + gP atr ) / 2 we restore approximately the isospin- 
selfconsistency of the total residual p-p interaction. 

• The p-p spin-isovector strength g^ p is chosen to reproduce the experimental value of M 2v . 
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Figure 1: The calculated gl -dependence of M 2u in both discretized and continuum QRPA. 



We perform the first calculations of M 2u and M 0u within the continuum- QRPA for 130 Te. Also we 
compare the obtained results with those calculated within the usual, discretized, version of the QRPA 
in order to to see the influence of the single-particle continuum. The chosen BCS basis contains 22 
levels (oscillator shells N = 1 -j- 5) that includes all bound s.p. states for neutrons and all bound s.p. 
states along with 6 quasistationary states for protons. The fitted values of the strength parameters f^ h 
and gp p are given in Table I for both discretized and continuum version of the QRPA. 

In Fig. 1 the calculated ^-dependence of M 2u is plotted. Note, that both (3~ and (3 + branches to 
construct the 2i/ / 3/?-amplitude are calculated for 130 Te, so we adopt here the same approximation as 
in [HI EI - 

The calculated values of M 0u are given in Table I for both versions of the QRPA and gA = 1.25. Note 
that the two-nucleon short-range correlations are included in the calculations in terms of the Jastrow 
function. At the same time, the higher order terms of the nucleon current are not considered (they 
usually reduce M 0u by about 30%, see, e.g. [2]). 

The contributions of the multipoles up to J 1 = 6" (L = 4- 5) are included in the calculations 
of M 0u . Note that the QRPA itself as a long- wave approximation is barely suitable to describe the 
multipole contributions with L > 5 (they contribute in total about 10% to M 0u ), they are completely 
dominated by the short-range behavior of the wave function. The partial multipole contributions to 
the calculated M 0u, s are given in Fig. 2. 

4 Conclusions 

In the article a continuum- QRPA approach to calculation of 2z//?/3- and Oz//5/3-amplitudes has been 
formulated. For 130 Te a regular suppression (about 20%) of the (L > 2)-multipole contributons to 
M 0u has been found which can be associated with additional ground state correlations appearing from 
the transitions to collective states in the continuum. At the same time the total M 0u for 130 Te gets 
suppressed by about 20% as compared to the result of the discretized QRPA. As the nearest perspective 
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Figure 2: Multipole decomposition of M 0u calculated within the both versions of the QRPA. 
J 71 " is the angular momentum and parity of the intermediate states. 



we are going to perform a systematic analysis of other double-decaying nuclei within the cQRPA. 
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